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Abstract The first generation of gravitational wave interferometric detectors 
has taken data at, or close to, their design sensitivity. This data has been 
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Fig. 1 Typical spectral density of calibrated noise for the LIGO interferometers and GEO 
during the S5 run and for Virgo interferometer during the VSR1 run. Also displayed are the 
design sensitivities for the 4 km LIGO interferometers and for Virgo. 



1 Introduction 

The first generation of gravitational wave (GW) interferometric detectors has 
reached an unprecedented sensitivity to GW signals. 

Six interferometric detectors completed commissioning activities and ac- 
quired scientific data over the last years. Three of them constitute the Laser 
Interferometric Gravitational- Wave Observatory (LIGO), a joint Caltech-MIT 
project supported by the National Science Foundation [Tj . They are situated 
in the USA, one in Livingston, Louisiana (LI) and two, which share the same 
facilities, in Hanford, Washington (HI, H2). LI and HI are interferometers 
with arms of length 4 km, whereas H2 is a 2 km interferometer. Their first 
Science Run began in September 2002. Since then, four other science runs 
took place; between November 2005 and September 2007 they operated at 
their design sensitivity in a continuous data-taking mode in the 5th science 
run (S5). 

The Virgo detector is a joint project of the CNRS and INFN, operated 
by the Virgo Collaboration at the European Gravitational Observatory - a 
CNRS-INFN joint venture with the mandate to build the interferometer [3]. 
Virgo is a 3 km interferometer near Pisa (Italy). Its first Science Run (VSR1) 
took place in May 2007-September 2007, with a sensitivity comparable to 
LIGO instruments above around 400 Hz and a better sensitivity below 30 Hz. 
(see Fig. 1}. 
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Two additional interferometric detectors are part of the network of GW 
observatories: the Japanese TAMA project built a 300m interferometer near 
Tokyo, Japan [3] and the German-British GEO project built a 600m interfer- 
ometer near Hannover, Germany [3]; early in its operation GEO joined with 
the LIGO Scientific Collaboration. The current sensitivity of TAMA and GEO 
is not comparable to LIGO and Virgo, and the latest scientific observational 
results are based on the analysis of S5 /VSR1 data. 

The LIGO S5 run collected a full year of triple detector coincidence inter- 
ferometer data, whereas Virgo ran in science mode for about 110 days during 
VSR1. One of the most promising candidates for a GW detection is thought 
to be the signal originated during the coalescence of two neutron stars. The 
distance at which a system of two neutron stars of 1.4 solar masses can be 
detected by a gravitational wave observatory, with a signal to noise ratio of 8, 
has thus become a standard measure of the detectors' sensitivity. During the 
S5/VSR1 run, HI reached a maximum sight distance of 35 Mpc, corresponding 
to a spherical observable volume of radius of 15 Mpc. 

LIGO and Virgo signed an agreement to jointly analyze data from the four 
LIGO- Virgo interferometers collected after May 2007. The Data Analysis ac- 
tivity is run by four LIGO- Virgo search physics groups with different scientific 
targets: Burst, Compact Binary Coalescences, Continuous Waves and Stochas- 
tic Background, each aiming at the detection and characterisation of different 
sources of GWs. Even if this division is somehow arbitrary, it has become the 
accepted standard in GW data analysis. 

Currently, a new joint run, S6/VSR2, is on-going. The goal is to improve 
the sensitivity of the detectors, Enhanced LIGO and Virgo+, by about a factor 
2 over the course of the run, resulting in an observable volume of the universe 
of about an order of magnitude larger. 

In the next sections of this article a general introduction to the methods 
applied and observational results obtained for the different targets will be 
presented. 

2 Bursts 

2.1 Introduction 

Gravitational wave bursts are taken to be short signals, duration < Is, with 
a large variety of possible waveforms, either wide or rather narrow-band in 
the frequency domain. GW bursts historically correspond to poorly modelled 
signals, due to the complex physics involved in the description of the sources, 
such as core collapse supernovae or mergers of two compact stars (black holes 
and/or neutron stars). Hence analysis methods seeking GW bursts must in 
general stay robust with respect to the expected waveforms. Moreover the 
benefit of running robust analysis methods is to stay wide open to the unex- 
pected, something which is fine when in a discovery mode. This constraint can 
be relaxed somewhat if some additional information about the source is avail- 
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able from other astrophysical observations. Also, spectacular recent progress 
in numerical relativity has allowed for accurate modelling of GW sources. For 
the most important burst sources we have now at least an idea of what the 
emitted waveform looks like. 

2.2 Sources of GW bursts 

Core collapse supernova — the collapse of a massive star and neutron star 
formation — has long been envisaged as one of the first GW sources. Modern 
models now include 3D full relativistic hydrodynamical simulations together 
with detailed micro-physics [3EEIE] . They predict the emission of GW signal 
during the fast gravitational collapse and bounce with a duration ~ 1 ms and 
peak amplitude h ~ 10 -21 for sources located at 10 kpc. A longer (hundreds 
of ms) and in some cases stronger signal is also expected after the bounce due 
to hydrodynamical instabilities around the proto-neutron star. 

Another source of GW bursts is the merger of two compact stars (black 
hole and/or neutron star) which generally results in a black hole, surrounded 
by a torus of matter if at least one of the progenitors is a neutron star. Progress 
in numerical relativity in recent years allows a clear picture of the merger of 
two black holes and a good prediction for the waveform [5UTU] . Good progress 
has been made too in modelling the merger of two neutron stars [11] , which 
is an even more complex problem due to, for example, the role of the neutron 
star equation of state. 

Core collapse supernova and merger of neutron stars are thought to be 
the origin of, respectively, long and short Gamma Ray Bursts (GRB). In both 
cases, a black hole will be formed and surrounded by a disk of matter. The 
accretion of matter onto the black hole leads to the formation of relativistic 
jets and the emission of gamma rays. What is interesting in this situation is 
that a GW signal and a GRB signal are emitted within a short delay. This 
may considerably enhance the detection capability, as explained below, and 
give deeper insight into the GRB mechanism. 

Soft Gamma Repeaters (SGR) sporadically emit brief and intense bursts 
of soft Gamma Rays and could also be a good source of GW bursts. SGRs 
are produced by highly magnetized neutron stars which undergo deformations 
(starquakes) that could excite the modes of the star and then emit gravita- 
tional waves at about the same time as the gamma emission. The expected 
waveform is typically a damped sinusoid ( "ringdown" ) but the parameters de- 
pend on the equation of state of the nuclear matter inside the star and are not 
accurately known. Like in the GRB case, the gamma ray emission permits a 
search for the GW bursts in a narrow coincidence time window. 

2.3 Detection methods 

Based on the classification of expected sources given in the previous section, 
there are basically two ways of searching for GW bursts: a general one making 
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very little or no assumption about the signal direction, waveform, time of 
arrival, etc., and a second where we take benefit of an external trigger (for 
example a GRB) for which source location in the sky and timing are known. 

The main steps of the analysis are trigger generation, efficiency estima- 
tion and accidental background estimation, whose principles are the same for 
both all-sky all-time unmodelled searches and triggered searches. Even with 
increasingly accurate predictions for GW burst waveforms, it is still mandatory 
that burst algorithms remain as robust as possible against the possible variety 
of waveform. Indeed, generic signals such as Gaussian pulses or Sine-Gaussian 
signals are used in parallel with astrophysical waveforms (collapse, mergers) to 
test the efficiency and robustness of trigger generators. The up-to-date trigger 
generator algorithms (or analysis "pipelines" ) are time- frequency or equivalent 
methods looking at local excess of energy (in time and frequency) in the cali- 
brated h(t) time series [T2lfT3l[T4U15lH6] . These methods can be fully coherent 
[T2lfT3. 16 taking full benefit of the existing network of detectors (LIGO and 
Virgo) or coincident with some coherent follow-up when looking for possible 
event candidates. A coherent method has the advantage that it can faithfully 
reconstruct the burst signal and the position of the source in the sky. Of course 
a coincident analysis can also reconstruct the sky location of some burst event 
candidate. In both cases accurate sky position can be reconstructed if the 
three LIGO and Virgo interferometers are in operation and in science mode 
at the same time. In the case of the current LIGO- Virgo network, the typical 
position angular accuracy is of the order a few degrees, depending on the GW 
direction with respect to the detectors plane [17llT8j . 

The efficiency estimation of the searches is performed by injecting signals 
in the calibrated time-shifted output data of the detectors. As already men- 
tioned, these signals can be astrophysical signals from core-collapse or merger 
simulations (see for example the search for the merger and ring-down phases 
of binary black hole coalescences with Virgo C7 commissioning run data [19] ) 
or generic signals. The so-called Sine-Gaussian signals (sinusoids with Gaus- 
sian envelope) are particularly interesting since they can span the detector 
bandwidth and give sensitivity estimates for the whole range of accessible fre- 
quencies. A number of such signals are then added to the calibrated strain of 
each detector after time sliding. For the all sky searches, they are generally 
distributed along the entire data in order to average the detectors' noise non- 
stationarities. For each injected signal the amplitude described by the h rssi 
which corresponds to the total GW signal energy: 



to the detector sensitivities, especially when the injected waveform is well 
localized in frequency. The final result for each signal is the efficiency curve, 
detection efficiency vs. /i rss . Usually the h rss values at 50% or 90% efficiency 
are chosen for setting the upper limits. 




(1) 



The 



h rss amplitude has units of Hz 1 / 2 , meaning we can compare it directly 
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Another important aspect of the burst analysis is the background estima- 
tion and the rejection of loud outliers (glitches) using data quality flags and 
event by event vetoes. The quality of the detectors' data is studied at the in- 
dividual detector level. Short duration noise transients can mimic a GW burst 
especially if their rate is large enough to induce a number of coincident events 
between the detectors. Data quality flags are used to indicate that one inter- 
ferometer is not working properly during some period, for some known reason. 
Event-by-cvcnt vetoes arc defined when an excess of coincidence is found be- 
tween the gravitational strain channel triggers and auxiliary, environmental 
channel events. The above triggers are generally not obtained with the search 
pipelines but by other algorithms which are required to be much faster since 
they run on hundreds of auxiliary channels. The procedure allows us to veto 
only short time intervals (sub-second), thereby saving the observation time 
of the detectors. Once the data is cleaned by use of data quality flags and 
vetoes it remains to estimate the background, i.e., the number of accidental 
coincidences. As the noise in each detector is independent, this is done by time 
shifting the detector data with respect to the others and processing all these 
time shifted data sets with the search algorithms. Relative time shifts are cho- 
sen to be much larger than the light times of flight between the detectors and 
much larger than the expected signal durations. 

There is a little difference for triggered searches. In this case, the data 
are split into an on-source region where the GW burst is expected and a 
background region where the noise is expected to be statistically similar to 
the one in the on-source region. The background region generally consists of 
a few hours from each side of the on-source region. All the efficiency and 
background studies are performed with the background data set. 

All the cuts of the searches are then defined with these shifted data sets and 
applied in a blind way to the original unshiftcd data set. Events passing all the 
cuts (if any) are then detection candidates and can be examined more closely. 
If no candidate is found the result of the search can be turn into exclusion 
plots, for instance rate of event vs. h Tm in the case of all-sky analysis. The 
frequentist upper limit for the rate of events with confidence level p is for 
example [20] 

i e(/l rSS ) 

where T is the observation time and e(h ISS ) is the detection efficiency for the 
considered waveform with amplitude h rss . 



2.4 Recent results 

The search for GW bursts has not yet yielded a detection and most published 
results give upper limits on strain (h rss ) or rate. When an astrophysical inter- 
pretation for a particular source is possible then these limits can be converted 
into meaningful astrophysical bounds. 
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The most recent all sky search for unmodelled GW bursts has been com- 
pleted using data from the hrst calendar year of the S5 LIGO run. The analysis 
has been split in two parts, a search for low frequency signals in a 64-2000 Hz 
band [51] and a high frequency search in a l-6kHz band [55] . The observation 
time T ~ 270 days sets an upper limit about 3.6 events/year for the rate of 
events at 90% confidence level (p = 0.9 in the formula above) in the lower 
frequency band. In the higher frequency band the network live-time is only 
T ~ 161 days, resulting in an upper limit for the rates of events about 5.4 
events/year at 90% confidence level. The strain limits depend on the details 
of the injected signals for efficiency estimates. The typical h rss limit for Sine- 
Gaussian signals and most Gaussian pulses is below 10 -21 Hz -1 / 2 in the low 
frequency band while it is up to a few 10 -20 Hz~ 1//2 in the high frequency 
band. Not surprisingly the search sensitivity follows the detector sensitivities 
which are better at low frequencies than at high frequencies. 

Another recent result concerns the search for a GW burst associated with 
GRB 070201 [23], detected by gamma-ray satellites. Those satellites found 
that the error box for the position of the GRB is centered at about 1 degree 
from the center of M31 (Andromeda) and overlaps the spiral arms of the 
galaxy. Andromeda is the closest spiral galaxy (at about 760 kpc) and an 
event possibly occurring at such a short distance would be outstanding. The 
GRB was a short one and likely progenitors for short GRBs are binary neutron 
star or neutron star-black hole coalescences or flares from SGR. The analysis 
for coalescence waveforms is presented in the next section. An unmodelled 
search for a GW burst in association with the GRB, yielded an upper limit for 
the radiated GW energy of about 4.4 x 10~ 4 M Q c 2 for GW bursts lasting less 
then 100 ms for isotropic emission occurring at the LIGO peak sensitivity near 
150Hz. This does not rule out the possibility of a SGR giant flare in Andromeda 
galaxy. Other searches for GW bursts associated to GRB have previously been 
pubblished by the LIGO collaboration [24] and Virgo collaboration 25J. 

Two searches for GW bursts possibly associated with SGRs have been 
recently published [551127] . The first one targeted the SGR 1806-20 and the 
SGR 1900+14, including a giant flare episode of SGR 1806-20 which occurred 
in December 2004 (during the LIGO "Astrowatch" period prior to the S4 
run) and a storm episode of SGR 1900+14 which occurred in March 2006 
(during the LIGO S5 run) . Using a set of different waveforms the search was 
able to set upper limits at 90% confidence level on the isotropic GW emitted 
energies in the range 3 x 10 45 to 9 x 10 52 ergs for a source located at 10 kpc. 
These upper limits depend on the detector sensitivities and antenna patterns 
at the time of the Gamma emission, on the loudest event in the on-source 
region and the injected waveforms for efficiency estimation. It is worth noting 
that some theoretical models predict maximal GW emission energy as high 
as 10 49 ergs. This is well in the range of the SGR analysis sensitivity. The 
second paper is a re-analysis of the SGR 1900+14 storm of March 2006 with a 
different ("stacking") method [27]. The gain in sensitivity is about one order 
of magnitude with respect to the first analysis [26] . 
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A more exotic analysis is the search for GW bursts emitted by cosmic string 
cusps. The first reported results used the 2005 LIGO data (S4 run) [28]. The 
search used matched filtering as the predicted signals have simple and well 
parametrized waveforms. Upper-limits have been set for the rate of events and 
for cosmic string parameters such as string tension loop size or reconnection 
probability. These limits are not competitive with the ones obtained by other 
cosmological observations like indirect bounds from Big Bang Nucleosynthesis, 
but analysis of the S5 data (much longer data taking and with sensitivity twice 
as good) might surpass current limits in a some portion of the cosmic string 
parameter space. 



3 Compact Binary Coalescence 

3.1 The Binary Coalescence Waveform 

Coalescing binaries comprised of black holes and/or neutron stars are ideal 
sources for gravitational wave detectors. During the latter stages of its evolu- 
tion, a binary emits gravitational radiation as the two component stars slowly 
spiral inwards before finally merging to form a single object which settles down 
to equilibrium. Indeed, the emission of gravitational waves from binary inspi- 
ral has been indirectly detected through observations of binary pulsars [29] . 
although the current gravitational wave frequency is too low to be observed 
in terrestrial gravitational wave detectors. As the inspiral progresses, how- 
ever, the frequency and amplitude of the gravitational waves increase. During 
the final seconds or minutes of inspiral and merger, the gravitational radiation 
emitted by systems with a mass between one and several hundred solar masses 
will lie in the sensitive band of ground based gravitational wave detectors. 

The precise form of the binary coalescence gravitational waveform depends 
sensitively upon the parameters of the binary, most notably the masses and 
spins of the binary components. The eccentricity of the orbit will affect the 
emitted waveform but, in most cases, it is expected that the binary will have 
circularized before entering the sensitive band of ground based detectors [30l 
[3Tj . Historically, the binary coalescence has been split into three parts: a 
slow inspiral, a highly relativistic merger and ringdown to a final equilib- 
rium state. Different techniques are used to calculate the waveform in each 
of these regimes. Depending upon the mass of the system, different stages of 
the evolution will emit gravitational waves at the sensitive frequency of the 
detector. 

When the components of the binary are widely separated, the orbit decays 
slowly due to energy emitted in gravitational radiation, and the waveform 
sweeps slowly upwards in both frequency and amplitude. During this inspiral 
phase, theoretical waveforms calculated within the post-Newtonian framework 
[32j are expected to provide an accurate representation of the gravitational 
waveform. The post-Newtonian waveforms derived to date are sufficient for 
binaries comprised of neutron stars or low mass black holes (up to about 
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IOMq), as the merger will occur above the most sensitive frequency band of 
the detector. 

It is expected that the end product of a binary system with total mass above 
2.0M© will be a single, perturbed black hole which rings down, by emission of 
gravitational radiation, to an equilibrium conguration. Since stationary black 
holes are fully characterized by their mass and angular momentum, the exci- 
tations of higher multipoles, and in partiular the quadrupole, will be radiated 
gravitationally |33j . The frequency and damping time for each of the ringdown 
modes depends upon the mass and angular momentum of the black hole and 
can be calculated analytically within the framework of black hole perturba- 
tion theory [54T 35 . Ringdown waveforms will lie in the sensitive band of the 
detector for black holes with mass greater than around 1OOM0. 

The dynamical merger of two black holes can only be modelled using full 
general relativistic calculations. Recent breakthroughs in numerical relativity 
have, for the first time, enabled the calculation of the gravitational wave- 
form emitted during merger [T0]l9"]. Currently, several groups are capable of 
numerically evolving two black holes through their final orbits, merger and 
ringdown [10, 36, 37J . The waveform for binaries whose components have com- 
parable mass, and are non-spinning, is well characterized |38 [l39lFTO ] . There 
have also been successes modelling the merger of neutron star binaries and 
neutron star black hole binaries [11,41,42, 43J. Numerical investigations of the 
full parameter space of compact binary coalescence waveforms are ongoing. 



3.2 Search Methods 

As described above, the gravitational waveform for coalescing binaries is well 
modelled analytically and numerically. Signal processing theory [44] advocates 
the use of matched filtering to extract known signals from Gaussian noise. 
Matched filtering provides the backbone of searches for coalescing binaries [HI 
[46] . However, two substantial challenges remain: searching over the large pa- 
rameter space of coalescing binary signals, and dealing with non-stationarities 
in the data. 

The full binary coalescence waveform depends upon as many as seventeen 
parameters and it remains a challenge to search efficiently over the full pa- 
rameter space. While some parameters, such as the amplitude and coalescence 
phase of the waveform, can be extracted using analytical techniques, others 
can only be searched by repeatedly evaluating the matched filter at numer- 
ous points across the parameter space. This is facilitated by creating a bank 
of template waveforms, subject to the condition that for any candidate sig- 
nal only a small fraction of the signal (typically 3%) is lost due to filtering 
with a mis-matched waveform [47,48,49j. This method works well for binaries 
with non-spinning components. However, the parameter space of spinning bi- 
naries is considerably larger. Several methods of attacking this problem have 
been proposed, including phenomenological waveforms |50p51j. restrictions to 
binaries with a single spin [52] and Markov Chain Monte Carlo [551153] tech- 
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niques. However, the increase in both computational cost and the background 
rate associated with covering the spin parameter space has, to date, rendered 
this unfeasible. Indeed, it was been shown that searching for spinning bina- 
ries with non-spinning waveform templates provides comparable sensitivity 
to the currently available spinning searches — the benefits of using the im- 
proved, spinning template model are negated by the increase in false alarms, 
particularly in real data |55j . Investigations of the use of spinning templates 
in gravitational wave searches continue. 

The data from gravitational wave interferometric detectors contains a sig- 
nificant number of non-stationary transients caused by various environmental 
and instrumental sources. These reduce the sensitivity of a matched filter 
search as loud noise transients will produce a large signal to noise ratio, even 
if they do not match well the gravitational wave binary coalescence signal. 
Many techniques have been developed to has mitigate the effect of these noise 
transients in the data. Among the most powerful are: data quality and veto 
tests which flag times of poor data quality or use auxiliary channels with 
known couplings to the gravitational wave channel to remove times of poor 
data [55] , as described in Section 12.31 coincidence tests which require that a 
signal be observed, with consistent parameters, at widely separated sites [56] : 
signal consistency tests which compare the observed signal in the detector to 
the predicted waveform [571158] : the use of improved ranking statistics which 
better separate the foreground and background by taking into account addi- 
tional information, over and above the signal to noise ratio of the candidate 
[60] . By making use of these additional tests, searches for binary coalescences 
are approaching the theoretically predicted sensitivity in Gaussian noise. 



3.3 Search Results and Future Prospects 

Gravitational wave data from the GEO, LIGO, TAMA and Virgo interfero- 
metric detectors have been analyzed for coalescing binary signals. To date, no 
gravitational wave signal has been observed. Consequently an ever improving 
set of upper limits has been placed on the rate of binary coalescence as a 
function of the mass of the binary. Upper limits have been derived for systems 
ranging from neutron star binaries through to intermediate mass black hole bi- 
naries. Here, we recap the latest results and compare them with astrophysical 
predictions. 

Astrophysical estimates of binary neutron star coalescence rates can be 
derived from observations of binary pulsars in the galaxy. These rates are 
extrapolated to the local universe under the assumption that the rate of binary 
coalescence follows the star formation rate in spiral galaxies, which is obtained 
from measurements of the blue light luminosity of galaxies [61]. Thus, results 
are quoted per L w per year, where 1L W = 10 10 times the solar blue luminosity 
and, for reference, the Milky Way is approximately 1.7Lio- The predicted rate 
for binary neutron star coalescence is 5 x 10 _5 L 10 1 yr _1 , although the rate 
could plausibly be as much as an order of magnitude larger [62] . To date, 
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there are no direct observations of black hole-neutron star or black hole-black 
hole binaries. Thus, rate estimates are based upon population synthesis, and 
yield realistic rates of 2 x 10~ 6 yr _1 Z^ for neutron star-black hole binaries 
[55] and 4 x 1CP 7 yr _1 Lj" for binary black holes [53] . In both cases, the rate 
could feasibly be as high as 6 x 1CP 5 yr -1 ^^ . 

The first eighteen months of the LIGO S5 data have been analyzed for grav- 
itation wave signals from coalescing binaries with a total mass less than 35M Q 
[65[l66j. Upper limits obtained by combining the result of this analysis with 
results from S3 and S4 provide the most stringent bounds on the coalescence 
rate from gravitational wave observations. The binary neutron star coalescence 
rate is restricted, at 90% confidence, to be less than 1.4 x 10~ 2 L ul 1 yr _1 [66] . 
This is a factor of thirty above optimistic rates, and several hundred above 
the best estimate of the rate. It is, however, interesting to note that the upper 
limit has improved by four orders of magnitude from the one obtained with 
LIGO's first science run [57]. For binary black holes, with component masses 
5 ± 1M , the 90% rate limit is 9 x 10 _4 ij~ 1 yr~ 1 , and for black hole-neutron 
star binaries, the limit is 4 x 10 _3 L ul 1 yr _1 [66]. These limits are again be- 
tween one and two orders of magnitude from the upper end of astrophysical 
predictions, and three orders of magnitude from best estimates. 

The most recent upper limits for binaries with a total mass greater than 
35-Mq were obtained using data from the fourth LIGO science run. For binaries 
of total mass between 30 and 80M Q a search with a phenomenological template 
family |68j of binary black hole waveforms gave an upper limit of ~ lL^yr -1 
[60] . For higher masses, a search for the ringdown portion of the signal yielded 
a rate limit for binary coalescences in the mass range 100 to 400Af Q of 1.6 X 
10 _3 L 1 Q L yr _1 . A search of the LIGO S5 data for black hole binaries with a total 
mass up to 100M Q is being pursued [69]. This search will, for the first time, 
make use of full inspiral-merger-ringdown coalescence waveforms obtained by 
enhancing the post-Newtonian inspiral waveforms with merger and ringdowns 
simulated numerically [38] . 

The coalescence of two neutron stars or a neutron star and a black hole 
is one of the preferred progenitor scenarios for short-duration GRBs [70]. By 
making use of the known time and sky location of observed GRBs, it is possible 
to perform a more sensitive search of the gravitational wave data. This has 
been done for GRB 070201, which was a short GRB localized in a region 
of the sky which overlapped the Andromeda galaxy [23]. The search yielded 
no evidence of gravitational waves, and allowed for the exclusion of a binary 
coalescence progenitor in M31 with 99% confidence. 

Recently data taking with the enhanced LIGO and Virgo+ detectors be- 
gan. It is hoped that these detectors will achieve a factor of two sensitivity 
improvement over the initial detectors, which translates to almost an order 
of magnitude increase in the volume of the universe that can be probed for 
binary coalescences. In the following years, advanced gravitational wave de- 
tectors will bring an order of magnitude increase in sensitivity over the initial 
configurations. At this stage, astrophysical estimates predict the observation 
of gravitational waves from tens of binary coalescences per year. 
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4 Continuous Waves 

4.1 Introduction 

Although neutron stars in coalescing binary systems represent a relatively 
well-understood population of putative gravitational wave sources, with well- 
defined gravitational luminosities and population statistics, the same neutron 
stars (and even isolated neutron stars) can in principle radiate a detectable 
amount of gravitational radiation well before coalescence. Radio and X-ray 
pulsar populations give us only a hint of the vast number (~ 10 9 — 10 10 ) of 
neutron stars that exists in the Galaxy. We currently see perhaps only one in 
a million neutron stars as a pulsar, but any neutron star, pulsar or not, can 
generate continuous quasi-sinusoidal gravitational waves through rotation. 



4.2 Sources 

Any non-axisymmetric spinning neutron star will generate gravitational ra- 
diation. Although the centrifugal deformation can be expected to make it 
significantly oblate (~ 1CP 4 ), this axisymmetric deformation will not itself 
generate gravitational radiation. Instead one requires the shape of the neu- 
tron star to be supported against relaxation to a fluid equipotential surface by 
a force, possibly an elastic stress force from the crust of the star, a magnetic 
force distorting the crustal shape or possibly a distortion caused by accretion 
or gravitational radiation-driven instabilities. A neutron star with its spin axis 
oriented towards us, at a distance d, with such a mass quadrupole moment Q 
around its axis of spin will generate a circularly polarised gravitational signal, 
at a frequency equal to twice the rotation rate v of the star, with amplitude 

16Gtt 2 v 2 n 

If the spin axis is inclined to the line-of-sight by an angle t, the radiation 
becomes elliptically polarised and the amplitude is reduced. It is often conve- 
nient to express Q as the product of an axial moment of inertia, I zz , and an 
effective equatorial ellipticity, 

e = — ~ Iyy , (4) 

Izz 

where I xx and I yy are the other two principal moments of inertia. There are 
clearly two very important questions here: i) what is the size of deformation 
that neutron stars could have, given their equation of state, crystalline struc- 
ture and physical environment, and ii) what deformations do they actually 
have. The answer to the first question depends critically on both the equa- 
tion of state of the neutron star (and whether it is indeed a neutron star or a 
quark star) and its crystalline structure [71072] . Recent work by Horowitz and 
Kadau |73j has indicated that the structure may indeed by highly crystalline, 



13 



with point defects rapidly squeezed out to give breaking stains of as much as 
0.1. This would allow self-supported equatorial ellipticities of perhaps 10~ 5 to 
1CT 6 . 

This second question is harder to address. Unlike in binary neutron star 
systems such as PSR B1913+16, PSR J0737-3039, PSR B1534+12 and PSR 
J1756— 2251, where the orbital evolution presents convincing evidence of the 
very early stages of a coalescence, we have no direct evidence of spin-gravitars 
(that is, neutron stars whose observed spin-down is well-modelled by gravita- 
tional braking) . Certainly in the case of equatorial deformation supported by 
crustal strength we cannot dismiss the notion that some neutron stars have 
perfectly annealed equipotential surfaces, with negligible axial quadrupole mo- 
ment and therefore essentially no gravitational luminosity. For example, the 
extremely low period derivatives of millisecond pulsars hint that these neutron 
stars at least show very little equatorial asymmetry. At some level all neutron 
stars will show deformation due to internal magnetic pressures, though this 
only becomes relevant for the strongest of magnetars [74]. However, the story 
is less clear for young and/or accreting neutron stars. For both these classes 
there are plausible mechanisms to supply both the energy and the deformation 
necessary for significant gravitational luminosity (see for example [7S176P77]). 

For an isolated neutron star we may postulate that the gravitational lu- 
minosity should be less than the rate of loss of rotational kinetic energy for 
a rigid body. In turn this defines an upper limit on the strain amplitude we 
could expect from rigid body gravitar with a spin-down rate i> as 

This "spin-down upper limit" falls below the 1-year strain sensitivity of both 
the initial LIGO and Virgo detectors for all but a small number of known ra- 
dio pulsars, making these pulsars unlikely first-detection candidates. However, 
only a tiny fraction of the neutron star population is seen as electromagnetic 
pulsars, leaving the possibility of a gravitationally luminous, but electromag- 
netically dim, population. 

4.3 Search methods 

There is an obvious sense in which it is easier to search for a continuous quasi- 
sinusoidal signal than a transient inspiral or burst signal. A long-lived signal 
can be re-observed (or even retrospectively observed using archived data) and 
confirmed as astrophysical. In addition, there are radio and X-ray pulsars that 
accurately trace the rotational evolution of around 100 pulsars whose gravi- 
tational signal would fall in terrestrial observing bands. However, the benefits 
of a continuous wave search stop there. An inevitable consequence of search- 
ing for a long-duration deterministic signal containing up to ~ 10 10 cycles 
is an exquisite sensitivity to its parameter values. Most obviously, a change 
in frequency corresponding to just one more or one fewer cycles during the 
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observation would represent an entirely different search, with a template for 
the expected signal that, as a matched filter, was insensitive to the original 
signal. It also becomes apparent that if the neutron star is spinning down, the 
number of spin-down templates necessary to cover the possible alternatives 
scales as the square of the observing time. Additionally, the doppler modu- 
lation of the received signal is sensitive to the position of the source on the 
sky. For year-long observations this angular sensitivity is approximately the 
gravitational diffraction limit of an aperture the diameter of the Earth's orbit 
about the Sun (~ 1 arcsec at typical frequencies). If we do not have the benefit 
of a radio trace of the neutron star's rotational evolution and sky location we 
are forced to perform a search over this parameter space, and it rapidly be- 
comes apparent that the parameter space is huge. Continuous-wave searches 
are by far the most computationally expensive searches that the gravitational 
wave community undertakes, and in its most general form the problem is (and 
always will be) fully limited by available computing power. 

No matter what its form, any coherent search method will improve its 
strain signal-to-noise ratio as 

snr cx {Sh/T)- 1 ' 2 , (6) 

where Su is the detector's (strain) power spectral density at the frequency of 
the signal and T is the observing time. However the overall sensitivity of a 
search is not solely dependent on signal-to-noise ratio. The more trials that 
are undertaken (i.e., templates that are searched) the greater the probability 
of random noise popping up to unluckily appear like a signal. The apparent 
signal-to-noise ratio indicative of a true signal is therefore several tens for 
searches that pick the strongest candidate from over a wide parameter space. 
One consequence of this is that any convincing signal must have a relatively 
high signal-to-noise ratio after only a relatively short coherent integration. 
This allows one to combine these short integrations incoherently (as powers, 
ignoring phase) without too great an impact in overall sensitivity and develop 
semi-coherent search methods which are computationally much cheaper |78j . 
The overall sensitivity does however only improve as the quarter power of the 
number of incoherently combined contributions. 

Ground-based CW search efforts have concentrated on variants of the 
above, from fully coherent long-timescale searches for gravitational wave sig- 
nals phase-locked to radio pulsars 79.80.81,82] to searches concentrating on 
non-pulsing targets [831184] and massively computational all-sky searches using 
a variety of semi-coherent techniques [85,86,87,88,89 . 

4.4 Search results 

Upper limits on the strength of continuous gravitational waves from both 
known and unknown galactic neutron stars have been made regularly since 
the first LIGO/GEO science run in 2002. As sensitivities and run lengths have 
improved, the limits have steadily dropped. The recent 23-month LIGO S5 run 
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had sufficient sensitivity to show that the Crab pulsar is not a gravitar (i.e., 
is not spinning down solely due to the emission of gravitational radiation). 
In itself this is no surprise - the overall energy budget of the Crab nebula 
and pulsar has to account for the nebula luminosity and expansion. However, 
the early S5 result (covering just the first 9 months of data) was sufficiently 
sensitive to show that less than about 3 percent of the spin-down luminosity 
of the Crab pulsar is due to gravitational emission [82] • The full S5 result, 
with a fully coherent search for gravitational emission from 116 known pulsars, 
including the Crab pulsar, is expected soon. Additional work is going on to use 
these and other newly-developed targeted algorithms to search Virgo VSR1 
data for emission from the Vela pulsar at ~ 22.5 Hz |90) . 

The S5 run has also resulted in the most sensitive "all sky" (i.e., survey) 
searches to date. Two early-S5 papers have already been published on this [HH1 
189] . The first comprised a semi-coherent search, incoherently adding 30-minutc 
demodulated power spectra (the "power flux method"). The second was an 
early result from Einstein@Home, a distributed screen saver application that 
currently attracts about 200000 users worldwide and returns ~ lOOTflops to 
search project. This method looked for coincident detections between multiple 
30-hour coherent searches. Both these all-sky searches returned strain upper 
limits of around 10~ 24 for a wide spectral band, and these are levels with 
real astrophysical significance. A simple argument, originally by Blanford but 
developed by Knipsel and Allen [91] indicates that for our Galaxy population 
and distribution of neutron stars, the loudest expected CW source would have 
a strain at Earth at about this level (under certain assumptions). Such a source 
would need to be within a few hundred parsecs of Earth. In addition to targeted 
and all-sky searches, more specialised "directed" search methods are also being 
used to tackle likely sky locations including globular clusters, the low-mass X- 
ray binary Sco-Xl, the galactic centre and supernova remnants. One such 
search, for gravitational emission from the X-ray point source at the centre of 
Cas A is nearing completion and has reported an expected sensitivity also in 
the range ~ 10~ 24 [92]. The Ligo and Virgo Collaborations have now developed 
a broad suite of algorithms and methods to tack a wide range of potential 
sources of continuous gravitational radiation, including all-sky searches for 
binary sources, and the full power of these will be applied to data form the 
current S6/VSR2 runs. 



4.5 Future prospects 

As with other searches that involve population statistics, the crude extrapo- 
lation holds that a factor r\ improvement in sensitivity will increase detection 
numbers by a factor ~ rf . Clearly the physical extent of the Galaxy places 
an upper limit on this, but that only becomes relevant for current all-sky 
searches when broadband sensitivities are a factor ^100 times their current 
values. Perhaps more important is a consideration of the types of neutron 
star that may be detectable in the future using instruments with an improved 
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low frequency response. Current detectors show good sensitivity only to rela- 
tively rapidly spinning pulsars, most of which are recycled millisecond pulsars 
with low observed spin-down rates and, probably, low gravitational luminos- 
ity. Young, glitchy pulsars are more common at gravitational frequencies below 
~ 100 Hz, with some of the most interesting, rapidly braked, sources closer to 
10 Hz, so the low-frequency wall is a particular challenge for future continuous 
wave gravitational observations. 



5 Stochastic Background 

A stochastic gravitational-wave background (SGWB) refers to a long-lived 
random GW signal. This is generally produced by a superposition of many 
unresolved sources, and can be characterized as cosmological or astrophysi- 
cal according to the epoch in which the GWs are generated. Cosmological 
backgrounds can be assumed to be approximately isotropic, unpolarized and 
stationary, while astrophysical backgrounds may have additional structure de- 
pending on the nature of their sources. 



5.1 Sources 

A convenient measure of the strength of a SGWB is the energy density in the 
GWs, per logarithmic frequency interval, in units of the critical energy density 
needed to close the universe!]] 

rtgw(/) = — #7 CO 
Pcritdln/ 

Cosmological models which produce a SGWB include amplification of quan- 
tum vacuum fluctuations during inflation [93,94,95], phase transitions (96] 
[97] . pre-big-bang models [98 l l99 lfT00] . and cosmic (super-)string models |101| 
I1021ll031ll04j . Standard inflationary models generate a backround of constant 
^gw(/) over many decades of frequencies, but the amplitude of such a back- 
ground is already bounded by cosmic microwave background observations to 
be f2 gw (f) < 10~ 14 |105j . Astrophysical GW backgrounds can be generated 
by unresolved superpositions of sources such as cosmic string cusps [104] . su- 
pernovae |106j . and neutron-star instabilities |107U108j . 

The most stringent indirect limit on a SGWB in the frequency range of 
ground-based detectors comes from a constraint on the total energy density 
present at the time of nucleosyntheis. This big-bang nucleosynthesis (BBN) 
bound limits the total energy density in gravitational waves to be 



yU g w(/) : - l<r"'(.V„-:J) (Si 



1 Note that p cr ;t = (3^qC 2 )/(87tG) depends on the value of the Hubble constant; it has 
become conventional to use the fiducial value Ho = 72km/s/Mpc when defining i?g W (/). 
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where the effective number N v of neutrino species at BBN is constrained to by 
(N v — 2>) < 1.44. [109] Note that this limit only applies to cosmological SGWBs, 
i.e., gravitational waves generated before the era of nucleosynthesis. 



5.2 Search Methods 

Since the amplitude of a SGWB will be much smaller than that of instrumental 
noise in a typical ground-based detector, one needs to exploit the expectation 
that while instrumental noise will be (predominantly) uncorrelated between 
independent detectors, the gravitational wave signals in a pair of detectors 
should have an average correlation 

(hiifTMf)) = \t(f J')7i 2 (/)5 gw (/) (9) 

where 712 (/) encodes the observing geometry (location and orientation of de- 
tectors 1 and 2, and in the case of an anisotropic background, the spatial 
distribution of the background) and S gw (f) is a one-sided power spectral den- 
sity for the SGWB which is given for an isotropic background by 

Sgw(/) = [(3 J ff 2 )/(10^ 2 )]/- 3 l2 gw (/) . (10) 

The standard search method |110j for an isotropic background cross-correlates 
the data from pairs of detectors using an optimal filter 

where Si^{f) are the noise power spectra for the two detectors and 6 gw (/) is 
the expected shape of the SGWB spectrum. The resulting search is sensitive 
to a background S gw (f) — Sn& gw (f) of strength 

/ 00 2\~ 1 / 2 

c<dctcctablc 

Sr ~[ 2T Jo df S 1 (f)S 2 (f) ) ■ (12) 

Note that the sensitivity of a cross-correlation search improves like the square 
root of the observing time T. Also, stochastic background measurements tend 
to be dominated by the low end of the available frequency range, because 
712 (/) oscillates with increasing / within an envelope whose leading term is 
oc / _1 and because Eq. [TO] means that a constant-.f2 gw (/) background has 

© S w(/) « r 3 . 

A cross-correlation search can also be used to search for an astrophysical 
background with a specified spatial distribution, e.g., a SGWB coming from 
one point on the sky [111] . More sophisticated techniques can be used to 
recover the spatial distribution of a measured background |112) . 
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5.3 Search Results 

The most stringent direct limit on f2 gw {f) was set using data from the S5 run 
of LIGO Livingston and LIGO Hanford [113] . which set the 95% confidence 
level upper limit of /2 gw (/) < 6.9 x 10~ 6 assuming J? gw (/) to be constant 
over the interval 41.5 Hz < / < 169.25 Hz. A SGWB of the excluded strength, 
confined to those frequencies, would contribute 9.7 x 10 -6 to the total value 
of fl. This limit is therefore more stringent than the BBN bound (Eq. [8]) and 
we have entered the era where ground-based GW detectors are placing new 
limits on gravitational wave backgrounds of cosmological origin. 

The previous limit from S4 data |114j . while less stringent by about an 
order of magnitude, already placed new restrictions on the parameters of some 
cosmic string models which generate GWs both before and after the era of 
nucleosynthesis. Additional searches of S4 LIGO data set limits on the strength 
of possible point-like backgrounds [115] and (by correlating LIGO Livingston 
data with data from the ALLEGRO bar detector) set a higher-frequency limit 
of f2 gw (915Hz) < 1.02 [116] . Correlation measurements using LIGO and Virgo 
data are expected to improve the high-frequency measurement [117] . Further 
searches for anisotropic backgrounds are also being conducted. 



6 Discussion 

The current search for GW covers multiple types of signals originating from 
different possible astrophysical events like core collapse of massive stars and 
neutron star formation, coalescing binary systems of neutron stars and black 
holes, non-axysimmetric spinning neutron stars and signals produced by a 
large collection of incoherent sources. The data acquired by the most sensible 
GW observatories, which are at present the LIGO and Virgo interferometers, 
are analysed applying different methods and strategies targeted to the identi- 
fication and characterisation of the signals emitted by these possible sources. 
Moreover, methods able to catch signals coming from unknown sources are 
also currently used. 

The analysis of the latest scientific data, acquired during the first part of 
the S5 /VSR1 run, did not provide any evidence of a possible detection. Upper 
limits on the rate of events and/or the strain amplitude h are then derived and 
compared to the astrophysical predictions. These limits are already valuable 
scientific results which reinforce and widen our knowledge of the astrophysical 
events involved. For example they are now approaching the plausible astro- 
physical values in the case of GW originating from binary coalescing systems 
and have already set a bound to the percentage of the spin-down luminosity of 
the Crab pulsar on the energy emitted in gravitational radiation. The energy 
density in a stochastic GW background around 100 Hz has been constrained to 
a limit which is more stringent than the the big-bang nucleosynthesis bound, 
the strongest indirect limit at those frequencies. From the analysis of the data 
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in coincidence with the GRB 070201 it has been possible to exclude the hy- 
pothesis of a binary merger in M31 as the progenitor of this event. 

Further results from the analysis of the full S5 / VSR1 run data are expected 
soon. Meanwhile, the analysis of the data from the S6/VSR2 run will proba- 
bly, if no GW detections are found, improve the current limits, owing to the 
possible improvement of sensitivity of LIGO and Virgo detectors. The sensi- 
tivity of the present and future GW detectors gives the possibility of studying 
astrophysical events jointly with other observatories, i.e. electromagnetic and 
neutrino observatories, likely bringing additional information on the physics 
of the sources and on their characteristics. The future searches will then open 
the possibility of performing a mature GW astronomy. 
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